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We investigate the interactions of closed strings in IIB matrix model. The basic interaction of 
tlie closed superstring is realized by the recombination of two intersecting strings. Such interaction 
is investigated in IIB matrix model via two dimensional noncommutative gauge theory in the IR 
limit. By estimating the probability of the recombination, we identify the string coupling in IIB 
matrix model. We confirm that our identification is consistent with matrix string theory. 
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I. INTRODUCTION 

IIB matrix model fl'| is considered as a candidate of 
the nonperturbative formulation of superstring theory. 
The relation between IIB matrix model and perturba- 
tive string theory is shown in Perturbative string 
theory is also contained in Dijkgraaf, Verlinde and Ver- 
linde's matrix string theory [3|. Strong coupling region 
of two dimensional supersymmetric Yang-Mills theory is 
described by the perturbative superstring theory. In the 
strong coupling limit, a free Green-Schwarz string theory 
is obtained. On the other hand, weak coupling region is 
described by the perturbative Yang-Mills theory. In ad- 
dition, there is an intermediate region which is described 
by type IIB supergravity solution in the large N limit ^] . 



The aim of this paper is to identify the string cou- 
pling gs in IIB matrix model. The hint for the identifi- 
cation comes from the matrix string theory. The gauge 
coupling of the two dimensional Yang-Mills theory has 
the \length]~^ of the worldsheet and it is related to the 
string coupling as gy^,/ = a' g^. The gs has the dimension 
[length] on the worldsheet. Thus, our task is to search 
for the dimensionful parameter on the worldsheet. 

Before searching for gs , we have to construct the world- 
sheets in IIB matrix model. They are constructed as two 
dimensional classical backgrounds [2,] in the IR limit. The 
string length is identified there and free multiple closed 
strings are obtained. Vertex operators of type IIA super- 
string arc constructed from IIB matrix model on these 
backgrounds in [3]. The relation between type IIA su- 
perstring theory and IIB matrix model is also proposed 
in [6] in a different way. 

In this paper, we consider the interaction of perturba- 
tive strings. The basic interactions are the transitions 
from two strings into one string or vice versa. These 
interactions are introduced in the formulation of super- 
string field theory in the light-cone gauge [7,]. 2 strings 
— > 1 string interactions are represented by the recombi- 
nation of intersecting strings locally. See figure 1. 




FIG. 1: Two strings intersect at an angle 
tion happens at the intersection point. 



The recombina- 
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Although this process locally represents the recombi- 
nation of 2 strings 2 strings, the final state is connected 
globally. Thus, after the recombination, we obtain a sin- 
gle closed string. In the matrix model, the instability 
of this system comes from the off-diagonal modes. The 
system with the larger intersection angle </> is more un- 
stable than that with the smaller intersection angle. In 
other words, the decay rate of two closed strings is de- 
cided by the intersection angle 0. We can replace the 
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angle (p with another parameter q which is related to (j) 
a.s (j) — 2tan~^q. Since we choose the horizontal axis 
in figure 1 as the coordinate of the worldsheet a and 
the vertical axis as the transverse scalar field 4>i{(j), the 
parameter q denotes the slope of the tilted strings. In 
the two dimensional gauge theory, q is the dimensionful 
parameter. We clarify the q dependence of the instabil- 
ity. Since the configuration decays into the single closed 
string through the recombination process, we estimate 
q dependence of the recombination probability. In the 
string perturbation theory, the recombination probabil- 
ity is proportional to gl to the leading order. 

In our investigation, we derive the action of multiple 
strings from IIB matrix model. We identify the unstable 
mode for the intersecting strings which are the solutions 
of this action. From this mode, we estimate the prob- 
ability of the recombination. By comparing the result 
to that of the perturbative string theory, we identify the 
string coupling in IIB matrix model. 

In section 2, we derive the action of multiple closed 
superstrings. In section 3.1, we carry out the fluctuation 
analysis for the intersecting closed strings. We identify 
the unstable mode. In section 3.2, we calculate the prob- 
ability of the recombination from the unstable mode. By 
comparing our result with that of the perturbative su- 
perstring theory, we identify the string coupling gs- In 
section 3.3, we estimate the probability of the recombina- 
tion in matrix string theory. Section 4 is devoted to the 
conclusion. In the appendix A, our notation of light-cone 
coordinates is given. In the appendix B.l, the differential 
equations for the fluctuation modes are solved. In the ap- 
pendix B.2, the Schrodinger equation which controls the 
time evolution of the probability is solved. 



II. THE EFFECTIVE ACTION OF MULTIPLE 
STRINGS 

The interaction among multiple strings is described in 
various ways in string theory. Perturbatively, transitions 
from two strings into one string or vice versa are the 
basic process. We aim to identify such interactions in 
IIB matrix model. Since this interaction is proportional 
to gs at the tree level, we can identify gs in IIB matrix 
model through this process. 

Let us start from the action 

S = -^Tr {^^[A>^,A^][A^,A,] + ]^4'T^[A^.^]^ ,(II.l) 

where tp is a, ten dimensional Majorana-Weyl spinor and 
A^{^ = 0, 1, • • • ,9) and -0 are NxN Hermitian matrices. 
By expanding the action pi.l|) around a two dimensional 
noncommutative background, 



we obtain an A/" = 8 two dimensional U{n) noncommu- 
tative gauge theory [8, 9, 10] 



S 



+2v;r'^i?^0 + 2v:;r,[0„^])^ , (11.3) 



where /i, = 0, 1 and i, j = 2, • • • ,9. 
The * product is described by 

a * & = exp ( a^„g ^ ) a{x + ^)bix + ry)|^=^=o , 



2 9^^577"- 



(11.4) 



where C^"^ is defined as the inverse of 6*^^. 

By taking the commutative (IR) limit, we obtain a 
commutative gauge theory 



5 = 



8vrg2 



+2^T''D~^^ + 2^r,[4, 7/-] + 0(9)) 



(II.5) 



Diagonal components are relevant degrees of freedom in 
the IR limit. We interpret the diagonal elements of the 
field 4>i in (jll.5|) as the coordinates of the fundamental 
strings. This two dimensional Yang-Mills theory is re- 
lated to a low energy effective theory of D-strings by the 
S-duality transformation, as is shown in fig. 2 of our 
previous paper [2]. 

We map the worldsheet coordinate from into x 
51 as 



z = xq + ixi 



r-\-it7 



(11.6) 



Since the only gauge invariant quantity is a set of the 
eigenvalues of the matrices (pi , if we go around the circle 
5^, the eigenvalues can be interchanged. See figure 2. 
We can consider a string of length n by the identification 
0i(o-) = (/)i(ct + 27m). 





FIG. 2: The string worldsheet coordinate. We map the coor- 
dinate from {xq,xi) into (t, cr). 



Since multiple free closed strings can be regarded as a 
moduli space of this theory, one can consider a configu- 
ration with multiple strings of various length in general 



{11.2) 



(/.,K) = (/>,(a'^ + 2™°) 



(II.7) 
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where n = J2a=i 

For our purpose, it is enough to analyze U{2) gauge 
theory since the recombination is a local problem which 
involves two strings. At the tree level, the amplitude of 
this interaction is proportional to Qs- After the coordi- 
nate transformation (|II.6p . the action (jll.Sp is mapped 
into 



S ^ -- 



87r.g2 



tr 2 



dTd(j\z\ 



+4 



[A+,A^] 



[(f>i,(f>j][(f>i,(l)j 



z 



[A^,ij]j+2^T,[^,,^]j{ll. 



Is- Since the recombination is the local interaction, the 
length of the string does not affect the interaction. Thus, 
for simplicity, we can consider the two closed strings with 
the equal length. We parameterize the integral region of 
(T as [0, 2ttw]. 

This action is valid when two strings coincide in the 
target space. If two strings intersect at a point which 
is indeed the situation we consider in this paper, then 
the off-diagonal elements of this action are relevant only 
near the intersection point. They are meaningful only 
two strings are very close to each other since otherwise 
these modes become massive. Thus, we can describe off- 
diagonal modes as local fields A^^ (r, a) which are valid 
near the intersection point with the small intersection 
angle. 

The definitions of light-cone coordinates d± and A± 
are summarized in the appendix lA II In the free string 
limit, we obtain light-cone Green-Schwarz superstring ac- 
tion which consists of marginal terms. In that case, every 
oscillator mode decomposes into a left mover and a right 
mover. 






FIG. 3: Two closed strings x and y recombine at some time 
and finally we obtain a single closed string z. 



After the field redefinition 
1/5" 




■4^R —i='4'R, 4>L 



and the rescaling : 





we obtain the effective action 



S 



ai -,12 



dr 



datr 
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+- 




--^\z\TpT,[(l),,il:] 



2^{T+D 

0{e)] . 



2 zz 



(11.9) 



(11.10) 



(11.11) 



The parameter ^ is identified with the string scale as 
a' = ^ in 0. 4>i has the target space dimension of 



III. RECOMBINATION 

A. Fluctuation analysis and the recombination 

In this section, by using the multiple closed superstring 
effective action (jll.lip obtained in the previous section, 
we analyze the closed string interactions. In the IR limit, 
diagonal components are relevant degrees of freedom. 
Thus, the theory becomes a free theory in this limit. 
The leading interaction in the perturbative string the- 
ory is the interaction between two closed strings, which 
is realized as the recombination. The recombination can 
be investigated by using the Yang-Mills theory in [ll| 
which is introduced as the low energy effective action of 
D-strings. 

Let us start from the two dimensional effective action 
pi.lip . The Bosonic part is written as 



S 



87rg2 J_ 



dr 



datr 



9 



2 zz 



OH' 



(111.1) 



In [3| , we identify the overall coefficient with g^^jy in 
the free string limit. Off-diagonal fields are regarded as 
local fields since we focus on the region where two strings 
are close to each other. 

The string coupling is very weak in the IR region. On 
the analogy of matrix string theory [3], the string cou- 
pling will behave gs (x We can interpret this relation 
as representing the equivalence between the IR limit and 
the weak coupling limit. 
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Supergravity solution of fundamental strings which is 
dual to our effective theory is described by [3l 



C/6 



4 077^4 



1 



6 o8^5i\r\-l/2 



-dU^ 



u 



2tt9 



2 df^^ , 



YM 



(III.2) 



We can read the scaling behavior of l/U^ from 

the metric. In the IR limit, the string coupling vanishes 
e"^ ~ t/"^ ~ l/x — > 0, which is consistent with our picture. 

It seems that this kind of running coupling behavior 
makes it difficult to treat the interaction. However, the 
recombination happens at a definite scale. For simplicity, 
we fix the scale of z as |z| = Zr = I and treat this process 
in the real time^ . We map the coordinates from z — e'^"'"*'^ 
into z = e'*^*"'"'^' by the analytic continuation r it. 
Then, the complex plane is mapped into the cylinder 
with the radius \z\ = 1. The Bosonic part of the action 
becomes 



+AD+(l>,D_i 



dt 



dati 



alp' 

f\ zz 



g- 



(III.4) 



This action is the two dimensional SU(2) Yang-Mills the- 
ory in the Lorentzian metric. Since we are considering 
the large winding number w, the worldsheet length of the 
string is very large. Thus, the solution of a single closed 
string might be represented locally as 



<P2 = y^a, 

^±=^ = 0, 0, = (i = 3,- • • ,9) , (III.5) 

in U{1) gauge theory. 

The solution which represents the intersecting strings 
is written as 




(i 3, - • • ,9) . (III.i 



For a generic value of Zr, the coupling constant is changed as 
3? ~* 9s/-^r • taking the following rescaling, 



we can absorb this factor. It is consistent with our claim (IIII.26I I 
9 ~ A- 



Two strings intersect at cr = 0. Each string is connected 
with each other at a point far from the origin a — Q, 
which cannot be seen in this local solution near the in- 
tersection point. Thus, if the recombination occurs at 
this point, we obtain the single closed string, g is a pa- 
rameter which is related to the intersection angle as 



tan ■ 



(III.7) 



q is the dimensionful quantity in the worldsheet sense 
and as we will see later, we relate it to the string coupling 
constant gs- If we take q <C 1, two strings are very close 
to each other in the target space, especially around the 
intersection point. Thus, our effective theory is valid. 

We consider the fluctuations around this background. 
We turn on the fluctuations a± and 



A+ = a+a"^ , A- = a-o^ , (III.8) 

since other fluctuations decouple from these fields at the 
quadratic level. f7'(i — 1, 2, 3) are the Pauli matrices. 
The quadratic lagrangian is obtained as 



—2{d+(p — iqa^a){d-(p — iqa^a) 



^\/2qa-(p + V2t 



qa+ip 



(III.9) 



Note that these fluctuation terms do not come from 
[4>i,4>j\[4>i,4>j\ terms. By choosing the gauge condition 



a 

V2 



(III.IO) 



the quadratic lagrangian is written by the parameter t 
and cr as 

+ (gacr)2 _ 2qaip)] (III.ll) 

In order to solve the equation of motion for the fluctua- 
tions, we expand the fluctuations by the mass eigenfunc- 
tions 



a{t,(T) = ^a(cr)C„(t) , 

n>0 



ri>0 



(III. 3) where C„(t) satisfy the equations 



(III.12) 



(III.13) 
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The differential equations are solved in the appendix 
IB 11 For the lowest mode n = 0, the eigenfunctions are 
calculated as 

Co{t) = Co-eMV^t) , 

5o(a) = ^o(0 -cxp(-|a2) , (III.14) 



with the eigenvalue 



(III.15) 



These eigenfunctions are localized near the intersection 
point (7 ~ 0. This mode is tachyonic and it causes the 
recombination 11]. 



B. The recombination probability 

Since the transition probability from two strings into 
the single string depends on the string coupling constant 
(7s, we can identify gs by our investigation. We estimate 
the recombination probability per unit time. Since we 
put two strings with the relative center of mass velocity 
w = as the initial condition, the recombination always 
occurs if we wait long enough. 

Before calculating the probability, we confirm that the 
lowest mode Cq (from now on, we denote it as C) indeed 
causes the recombination. During and after the recombi- 
nation, we investigate the time evolution of the tachyonic 
mode. Small fluctuations trigger the recombination and 
if we wait long enough, the recombination is over. 

The recombination can be seen by the diagonalization 
of the background with the off-diagonal fluctuations 



/ ^iqa)^ + C^it)e-i-' 

(III.16) 

Note that although this phenomenon locally represents 
the recombination of 2 strings 2 strings, the final state 
is connected globally. Thus, this recombination repre- 
sents the transition from two closed strings into the sin- 
gle closed string. We need to clarify the validity of the 
approximation since this geometrical picture comes from 
the fluctuation analysis. Our analysis is valid if the terms 
quadratic in the fluctuations are much smaller than the 
terms in the background. This condition is 

{qaf > C^e-^'^" . (III.17) 

Since qa^ ~ 0(1) as can be seen in (|III.19[) . we obtain 

(J > . (III.18) 



Our analysis is valid under this condition. 

The action quadratic in the C mode is obtained as 



S 



1 

8^ 



dt 

— OO <y —7TW 



da{{dtC{t)f + qC^{t)) 
xexp(-V) 

-^\jYqj ((atC(i))^ + qC\t)) . (III.19) 

After integrating the a direction, this action can be re- 
garded as a quantum mechanics of a particle moving in 
the inverse harmonic oscillator [l3|. As seen in piI.16|l . 
C{t) is related to the separation of recombined strings. If 
C takes a large value, then we can regard it as the signal 
for recombination. Once C obtains the large value, the 
recombination is over. The separation of strings grows 
and finally we obtain a single closed string. In the in- 
verse harmonic oscillator potential, the wave functions of 
the particle will keep spreading to the larger |C|. Thus, 
the recombination probability approaches 1 at t ^ cxd. 

The parameters in (jlll.lQP can be interpreted in terms 
of the quantum mechanics of the particle moving in the 
inverse harmonic potential as 



m 




(III.20) 

where m is the mass of the particle and lo is the frequency. 
The Schrodinger equation is given by 



.dip 



1 d'^ip 



(III.21) 



dt 4m 9C2 

As an initial condition at t = 0, we consider the wave 
function to be a gaussian which is labeled by the param- 
eter (j). The derivation of the solution of the equation 
pil.211) is discussed in the appendix IB 21 For large t, the 
wave function behaves 



V'(C,i) - (2/7r) 1/^6-1/2 exp 



1 



X exp 



0^ 



.(III.22) 



As discussed previously, once recombination happens, we 
obtain the single closed string as a final state. Since 
the fluctuation analysis is valid in the region ^/q 3> C, 
the geometric picture is also reliable in this region. We 
judge that the recombination has happened if the value 
of C(t) grows beyond ^ — uj. Thus, the recombination 
probability at a time t is estimated as 

/•OO 

P{i) = 2 / dC|7^(C,t)p 



= 1 - Erf(%/26-ie-*"w) 



1 - Erf [y^muj^ sin 2(j)e'*'^^ . (III.23) 
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One can confirm that at i ^ oo, P ^ 1. The recombi- 
nation probabiUty per unit time is calculated as 



dP 

dt ^jTi 

21/4 



— = -^^4^71^^ sin20exp (-4ma;^ sin20e - tuj] 



7r3/4 
X exp 



qy^ sin 2(f> 

^ sin2,^e-2\/9*-y^t^ . (III.24) 



2V27r 

In the small q limit, the probability is proportional to 
dP 



dt 



cx q 



(III.25) 



In the perturbative string calculation, this probability is 
proportional to g^. Thus, we identify 



(III.26) 



The higher order corrections come from the expansion of 
the exponential function. By expanding this function at 
a time t ~ 0{-^), and in the small q limit, we obtain 



(III.27) 



n=l 



which is consistent with the higher order corrections of 
the perturbative string if we identify q with 



C. The recombination in matrix string theory 

In this section, we estimate the recombination prob- 
ability of intersecting strings from matrix string theory. 
The action of matrix string theory is written by the two 
dimensional Yang-Mills theory 



-1 foo /'2'7rw 

S = ^ dr dai, {{lWs)F: 



2 2\ 1^2 



+2{D+4>,D+(t>,. + D-(t>,D^cp,) 

'■s9s 



(III.28) 



where we consider long intersecting strings with the large 
winding number w. By the identification = p-, the ac- 
tion (|III.28P is very close to the action pil.4p . The differ- 
ent point is the gs dependence which appears explicitly 
in pil.28[) . Thus, we perform the fluctuation analysis 
for the action pil.28[) and investigate the gs dependence. 
Since the action is the same apart from the gs depen- 
dence, if we consider the same solution, we obtain the 
same tachyon mode as in section 3.1. 



By regarding the tachyon effective action as the quan- 
tum mechanics of the particle moving in the inverse har- 
monic oscillator, the parameters gs and q appear in the 
mass and frequency of the particle as 



m = 




(III. 29) 



The only different point with respect to the analysis in 
section 3.1 is the gs dependence. The recombination 
probability at a time t is estimated from the previous 
calculation pil.24[) as 

P{t) = 1 - Erf 1^^452 y|gsin2^e-v^*j . (III.30) 

Thus, the recombination probability per unit time is 
dP 2I/4 



dt 7r3/4 
X exp 

By the rescaling 

we obtain 
dP 2I/4 



gsq\/sm2(t> 
1 



2V27r 



.gf q sin 20e^2tv^ - y/qt] . (III.31) 



t gst 



(III.32) 



dt 7r3/4 
X exp 



.9s 9 



■\/ sin 2(f) 



2V2tt 



sm2(l)e-^3st^ ^ .(III.33) 



This result is equivalent with the previous result pil.24p 
if we replace q qg'^. Thus, the identification we have 
derived in the section IIII Bl is consistent with matrix 
string theory. 



At a large time t ^ 0{ 



in the small q limit, the 



leading contribution is proportional to g^q. The higher 
order corrections which come from the expansion of the 
exponential part are 



(leading contribution) 



(III.34) 



which is consistent with the perturbative string theory. 

In [1^, the classical BPS solutions that interpolate 
between the initial and the final string configurations 
are constructed in matrix string theory. They interpret 
the amplitude of matrix string theory as the transition 
amplitude between initial and final configurations, and 
show that the leading contribution is proportional to 
gj^ where x is the Euler characteristic of the interpo- 
lating Riemann surface. Thus, they have reproduced the 
perturbative string amplitude from matrix string theory, 
which is consistent with our result. 
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IV. CONCLUSION 



A± are related to Aq and Ai as 



We have identified the string coupUng in IIB matrix 
model. We have constructed the classical solution of the 
strings in the action which is obtained from IIB matrix 
model. Starting from the configuration with the intersec- 
tion angle (j) = 2tan~^ q and no transverse distance and 
relative velocity, the recombination happens. This is trig- 
gered by the tachyonic fluctuations around the classical 
solution. After the recombination, we obtain the single 
closed string. We have estimated the probability of the 
recombination of two strings per unit time in this situa- 
tion. The recombination probability is also calculated by 
the perturbative string theory. The leading contribution 
is proportional to g'^. The higher order corrections are 
seen as J2'^=2 ^nds^- Comparing our results with these 
behaviors, we have identified the string coupling q ^ g-^. 

We have also estimated the recombination probability 
per unit time in the matrix string theory. The result 
shows the consistent behavior with that of the pertur- 
bative string theory. In the matrix string theory, the 
Yang-Mills coupling has the dimension [length]^^ in the 
worldsheet and it is inversely proportional to the string 
coupling gg. Thus, gs has the dimension [length] in 
the worldsheet. Since q has the worldsheet dimension 
[length]~'^ , the dimensionless parameter is qg"^. Since the 
probability is the dimensionless quantity, q should appear 
with g^ which is consistent with the result obtained from 
IIB matrix model. 

The parameter q represents the ratio between the tar- 
get space coordinate and the worldsheet coordinate 
a. Another parameter which has this kind of property is 
the velocity v. Although the recombination probability 
depends on the relative velocity u, we have put u = in 
this paper. It might be interesting to calculate the prob- 
ability of the recombination of intersecting strings with 
nonzero relative velocity. 
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APPENDIX A: NOTATION 



1. Light-cone coordinates 



d± are the derivatives with respect to 



1 



y± = ' 



and they are related to do and di as 

9o = — (— + —), di = —{^~- 

^ z z ' z z 



(A.l) 



(A.2) 



Ao^^{A+ + A^), Ai = -J=(^+-A_) . (A.3) 

The definition of F^g in (|IlTT|) is 

F,, = zd+{z-'A^)-zd-{z-^A+)-[A+,A^] 

= d+A^ - d-A+ - [A+,A^] . (A.4) 



2. Pauli matrices 



The Pauli matrices in piI.Sp satisfy the following rela- 
tion 



a\a^ = ie''''a'' , trier 



i\2 



1 



(A.5) 



APPENDIX B: CALCULATION 
1. The equation of motion for (|III.lH) 



The derivation of the eigenfunctions piI.14p is summa- 
rized in this appendix. 
By using the relation 

d+ + d- = -V2idt , 

d+-d- = -V2ida , (B.l) 

and imposing the gauge fixing condition 0+ = — a_ = 
we obtain the lagrangian pil.lip . 
The equation of motion for the fluctuation lagrangian 
(HIlTTI) is 



where 



O = 



o{ "S^;^} 1=0 



-{qcrY - d'^ -qada + q 
qada + 2q -92 + 92 



(B.2) 



(B.3) 



For the mass eigenvalues to^ which satisfy the free field 
equation 



the equation of motion is given by 

-{qa)'^+nil -qad^ + q \ ( a{a) 



qada + 2q + ml 



^(cr) 



(B.4) 



. (B.5) 



This differential equation is solved with the mass eigen- 
value 



ml = (2n- l)q 



(B.6) 
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For the lowest mode n = 0, the eigenfunctions are calcu- 
lated as 

Co(r) = Co • exp(7gi) , 

do{a) = (^o(o-) = exp (^-|cr^^ . (B.7) 

For general n, the eigenfunctions are 



We describe the momentum conjugate to the C as 11 

(B.IO) 



n = ^ = 2mC 
dC 



Hamiltonian is given by 



H = nc-L 



o„(cr) = -e" 



" n(n-2)...(n-j + 2) 

^ ^ 7'! 2n - 1 



4to 



(B.ll) 



x(j-l) U 



2 

(^nM = e""^ ^ (-1) 



. 4* n(n - 2) • • • (n - i + 2) 



2n~\ 



x(2n-j-l) (^a^l 
for n = 0, 2, • • • , and 



a„(cr) = -e '2 ^ (_i 



(B. 



, 0-1) / . , 

4 2 / J — 1 



X (n - 3) • • • (n - j + 2) cr 



2 

(^„(cr) = e""^ ^ (-1) 



4^— / J + 1 
n - 



x(n-3)..-(n-j + 2)(aJ^) (B.9) 



for n = 3, 5, • • • . 



2. Solving the Schrodinger equation HIII.2H) 



We summarize the derivation of the solution (jIII.22[) 
in this appendix. 



The Schrodinger equation (jIII.21|l is solved by the wave 
function of the form 

V'(C, t) = A(t) exp{-B{t)C^) , (B.12) 
if A and B satisfy the equation 



iA = —AB , 
2™ 

-B = OJ^ + ^B^ 
m 



(B.13) 



These are solved in llQl as 



B =^ JTLLU tan((/) — iujt) 

sin 2(j) — i sinh 2wt 

— muj , 

cos 20 + cosh 2ijjt 



(B.14) 



where 



,2 _ 



2mLu 

b = a(sin2</>)~i/2 _ 

We put the initial condition of the wave function as 
= 0) = (2^)-i/4(6cos0)-i/2 

X exp(— rnijjtan^C^) . (B.16) 

The parameter (j) controls the initial condition. 
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